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^ ■ A microscopic few-body description of near-threshold coherent photoproduc- 

tion of the r] meson on tritium and ^He targets is given. The photoproduction 
cross-section is calculated using the Finite Rank Approximation (FRA) of the 
nuclear Hamiltonian. The results indicate a strong final state interaction of 
the T] meson with the residual nucleus. Sensitivity of the results to the choice 
of the rjN T-matrix is investigated. The importance of obeying the two-body 
■ unitarity condition in the t]N system is demonstrated. 

^ : PACS numbers: 25.80.-e, 21.45.+V, 25.10.+S 

O 

O ■ I. INTRODUCTION 

Investigations of the ?7-nucleus interaction are motivated by various reasons. Some of 
^ , them, such as the possibility of forming quasi-bound states or resonances in the r^-nucleus 
pH I system, are purely of nuclear nature. The others are related to the study of the properties 
and structure of the 5*11(1535) resonance which is strongly coupled to the rjN channel. 

For example, it is interesting to investigate the behavior of the r^-meson in nuclear media 
where, after colliding with the nucleons, it readily forms the Su resonance. The interaction 
of this resonance with the surrounding nucleons can be described in different ways 0, 
depending on whether the structure of this resonance is defined in terms of some quark 
configurations or by the coupling of meson-baryon channels, as suggested in Ref. 0,^]. The 
estimation by Tiwari et al. |]^ shows, that in case of pseudoscalar rjNN coupling there is an 
essential density dependent reduction of the 77-meson mass and of the t] — t]' mixing angle. 

The importance of the influence of the nuclear medium on the mesons passing through 
it, was recently emphasized by Drechsel et al. [^. If this influence is described in terms 
of self-energies and effective masses, then in the process of vr-meson passing through the 
nucleus, "saturation" of the isobar propagator (or self-energy) takes place. This phenomenon 
manifests itself even in light nuclei @]. Similar ideas were discussed also in Ref. 0. In other 
words, the propagation of r^-mesons inside the nucleus is a new challenge for theorists. 

Another interesting issue related to the r^-nucleus interaction, is the study of charge 
symmetry breaking, which may partly be attributed to the rj — n^ mixing (see, for example, 
Refs. [pHTTH). In principle, one can extract the value of the mixing angle from experiments 
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involving ?7-nucleus interaction and compare the results with the predictions of quark models. 
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However, to do such an extraction, one has to make an extrapolation of the 77-nucleus 
scattering amphtude into the area of unphysical energies below the 77-nucleus threshold. 
This is a highly model dependent procedure requiring reliable treatment of the 77-nucleus 
dynamics. 

In this respect, few-body systems such as rjd, 77 ^He, and rj^He, have obvious advantages 
since they can be treated using rigorous Faddeev-type equations. To the best of our knowl- 
edge, so far only the simplest of these systems, namely the ri{2N) system, was considered 
within the exact AGS theory ]16[. 

A solution of the few-body equations presupposes the knowledge of the corresponding 
two-body T-matrices t^Tv and oS the energy shell. Due to the fact that at low ener- 
gies the 7] meson interacts with a nucleon mainly via the formation of the S'n-resonance, 
the inclusion of the higher partial waves {i > 0) is unnecessary. Furthermore, since the 
T]N interaction is poorly known, the effect of the fine tuned details of the "realistic" A^A^ 
potentials would be far beyond the level of the overall accuracy of the rjA theory. Indeed, 
in contrast to the well-established A^A^ forces, the t]N interaction is constructed using very 
limited information available, namely, the rjN scattering length and the parameters of the 
5'ii-resonance. Furthermore, only the resonance parameters are known more or less accu- 
rately while the scattering length (which is complex) is determined with large uncertainties. 
Moreover, practically nothing is known about the off-shell behavior of the rjN amplitude. 
It is simply assumed that all mesons should have somewhat similar properties and therefore 
the off-shell behavior of this amplitude could be approximated (like in the case of vr mesons) 
by appropriate Yamaguchi form-factors (see, for example, Refs. [p!2| - [T5|JT7|jr8[| ) . However, 
if the available data are used to construct a potential via, for example, Fiedeldey's inverse 



scattering procedure [|T9|, the resulting form factor of the separable potential is not that 
simple. The problem becomes even more complicated due to the multichannel character of 
the riN interaction with the additional off-shell uncertainties stemming from the vr-meson 
channel. 

In such a situation, it is desirable to narrow as much as possible the uncertainty intervals 
for the parameters of riN interaction. This could be done by demanding consistency of 
theoretical predictions based on these parameters, with existing experimental data for two-, 
three-, and four-body //-nucleus processes. This is one of the objectives of the present work. 
To do this, we calculate the cross sections of coherent 77-photoproduction on ^He and ^H 
nuclei and study their sensitivity to the parameters of rjN amplitude. 



II. FORMALISM 

We start by assuming that the Compton scattering on a nucleon, 

7 + A^ Ar + 7 , 

as well as the processes of multiple re-appearing of the photon in the intermediate states, 

-f + N^N + T]^'y + N^N + r]^... , 

give a negligible contribution to the coherent r^-photoproduction on a nucleus A. Then the 
process 
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+ + (1) 

can be formally described in two steps: at the first step, the photon produces the rj meson 
on one of the nucleons, 

-f + N + , (2) 

and at the second step (final state interaction), the rj meson is elastically scattered off the 
nucleus, 

1] + A-^ A + r] . (3) 

An adequate treatment of the scattering step is, of course, the most difficult and crucial 
part of the theory. Among the few-body systems r]d, rj^B., rj^He, and rj^He, so far only 
the simplest three-body one (rid) was considered in the framework of exact Faddeev-type 
AGS equations. The first microscopic calculations concerning the low-energy scattering of rj- 
meson from ^H, ^He, and ^He nuclei were done in Refs. [p0|-p6| where the few-body dynamics 



of these systems were treated by employing the Finite-Rank Approximation (FRA) of 
the nuclear Hamiltonian. This approximation consists in neglecting the continuous spectrum 
in the spectral expansion 

HA = Yj^n\i'n){i'n\+ coutinuum 

n 

of the Hamiltonian Ha describing the nucleus. Since the three- and four-body nuclei have 
only one bound state, FRA reduces to 

Ha-So\^o){H- (4) 

Physically, this means that we exclude the virtual excitations of the nucleus during its 
interaction with the t] meson. It is clear that the stronger the nucleus is bound, the smaller 
is the contribution from such processes to the elastic rjA scattering. By comparing with the 
results of the exact AGS calculations, it was shown p8[ that even for rjd scattering, having 
the weakest nuclear binding, the FRA method works reasonably well, which implies that 
we obtain sufficiently accurate results by applying this method to the rj ^H, rj ^He, and even 
more so to the t] ^He scattering. 

In essence, the FRA method can be described as follows (for details see Ref. [^). Let 



H = ho + V + Ha 

be the total rjA Hamiltonian, where ho describes free ?7-nucleus motion and 



A 



the sum of the two-body r^-nucleon potentials. The Lippmann-Schwinger equation 

A A 

T{z) = E + E ^^(^ - ^0 - Ha)-'T{z) (5) 

i=l 1=1 
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for the 77-nucleus T-matrix can be rewritten as 



T{z) = W{z) + W{z)M{z)T{z) , (6) 

where 

Miz) = Go{z)HaGa{z) , (7) 

Goiz) = (z - ho)'' , (8) 

GAiz) = iz-ho-HA)-' , (9) 

and the auxihary operator W{z) is spht into A components of Faddeev-type, 

W{z) = j:W,iz) , (10) 
1=1 

satisfying the following system of equations 

A 

W,{z) = U{z) + U{z)Goiz) J2 Wj{z) (11) 

with tj being the two-body T-matrix describing the interaction of the ?7-meson with the i-th 
nucleon, i.e. 

U{z) = V, + ViGoiz)U{z) . (12) 

It should be emphasized that up to this point no approximation has been used yet and 
therefore the set of equations (^fj^) equivalent to the initial equation (^. However, to 
solve Eq. we have to resort to the approximation (^) which simplifies its kernel (|^) to 

M(z) « , ff'X-^'l .. . (13) 

{z - ho){z -to- ho) 

With this approximation, the sandwiching of Eq. (||) between {ipo\ and l^po) and the partial 
wave decomposition give a one-dimensional integral equation for the amplitude of the process 
(0). 

The question then arises on how can a photon be included into this formalism in order 
to describe the photoproduction process (||). This can be achieved by following the same 
procedure as in Ref. [|29| where the reaction with A = 2 was treated within the framework 



of the exact AGS equations, and the photon was introduced by considering the r]N and jN 
states as two different channels of the same system. This implies that the operators ti 
should be replaced by 2 x 2 matrices. It is clear, that such replacements of the kernels of 
the integral equations ( [Tl| ) and subsequently of the integral equation lead to solutions 
having a similar matrix form 

Here t]'^ describes the Compton scattering, t^'^ the photoproduction process, and t^^ the 
elastic 1] scattering on the i-th nucleon. What is finally needed is the cross section 
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da 2 E^mA 



dn 9(27r)2 kr, + niA 



{kr„tPo\T'^\So + E,)\iJo,k,) (15) 
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of the reaction (|I|), where k^ and /c^ are the momenta of the photon and rj meson, E^ is the 
energy of the photon, the mass of the nucleus, and /i^a the reduced mass of the meson 
and the nucleus. 

However, it is technically more convenient to consider the ry-photoabsorption, i.e. the 
inverse reaction. Then the photoproduction cross section can be obtained by applying the 
detailed balance principle. The reason for this is that all the processes in which the photon 
appears more than once, i.e. the terms of the integral equations of type W^^MT^"^ or 
W'^^MT'^^ involving more than one electromagnetic vertex, can be neglected. Omission of 
these terms in results in decoupling the elastic scattering equation 

^ yyvv ^ w'^'^MT'^'^ (16) 

from the equation for the photoabsorption 

T'^'' = W^"^ + W^'^MT'^'^ . (17) 

Once the T*'^ is calculated, the photoabsorption T-matrix (|l^ can be obtained by integra- 
tion. 

Therefore, the procedure of calculating the photoproduction cross section ( 1T5|) consists 
of the following steps: 

• Solving the system of equations 

wr = tT + tTG^j2^T (18) 

for the auxiliary elastic-scattering operators W^^ . 

• Calculating (by integration) the auxiliary matrices W^^ from 

wr = tT+trG,jzw]^ . (19) 

• Solving the integral equation 

A A 

T^" = Y^W^'' + WJ'^MT'^'^ (20) 

i=l 1=1 

for the elastic-scattering T-matrix. 

• Calculating (by integration) the photoabsorption T-matrix 

A A 

T^^ = Y^W^'^ + Y^ WrMT"^"^ . (21) 

i=l i=l 



Substituting this T-matrix into Eq. (15) to obtain the differential cross section for the 
photoproduction. This is possible because the absolute values of the photoproduction 
and photoabsorption T-matrices coincide. 
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III. TWO-BODY INTERACTIONS 



To implement the calculation steps described in the previous section, we need the two- 
body T-matrices f'' and f^^ for the elastic ?7A^ scattering and the photoabsorption N{r], 'y)N 
on a single nucleon respectively. Furthermore, all equations (p^|2lD have to be sandwiched 
between {ipo\ and \ipo) (ground state wave function of the nucleus). Since at low energies 
both the elastic scattering and photoproduction of t] meson on a nucleon proceed mainly 
via formation of the Sn resonance, we may retain only the S-waves in the partial wave 
expansions of the corresponding two-body T-matrices. 



A. Elastic rjN scattering 

The problem of constructing an rjN potential or directly the corresponding T-matrix 
t"^^ has no unique solution since the only experimental information available consists of the 
^ii-resonance pole position Eq — iT/2 and the rjN scattering length a^jy. In the present 
work, we use three different versions of t^^ . 



1. Version I 

With the absence of any scattering data it is practically impossible to construct a reliable 
rjN potential. However, we can make use of the fact that the 5*11 -resonance is the dominant 
feature of the rjN interaction in the low-energy region, where the elastic scattering can be 
viewed as the process of formation and subsequent decay of this resonance, i.e. 

r] + N^Sn^N + r]. (22) 

This implies that in this region the corresponding Breit-Wigner formula could be a good 
approximation for the rjN cross section. Therefore, we may adopt the following ansatz 

t^^{k',k;z)=gik')T{z)gik) (23) 

where the propagator t{z) describing the intermediate state of the process (|22|) , is assumed 
to have a simple Breit-Wigner form 

= ^A-^ , (24) 

z — Eq + iT/2 

which guaranties that the T-matrix (^) has a pole at the proper place. The vertex function 
g{k) for the processes rjN ^ Sn is chosen to be 

g{k) = {e + aY' (25) 

which in configuration space is of Yukawa-type. The range parameter a = 3.316 fm~^ was 
determined in Ref. while the parameters of the S'n-resonance 

Eo = 1535 MeV - (mjy + m^) , T = 150 MeV 
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are taken from Ref. The strength parameter A is chosen to reproduce the ?7-nucleon 

scattering length a^Tv, 

A = ZTT ttriN ■ (26) 

the imaginary part of which accounts for the flux losses into the vrA^ channel. Here /x^at is 
the 'qN reduced mass. 

The two-body scattering length a^iv is not accurately known. Different analyses ||3^ 
provided values for o^at in the range 

0.27 fm < Re a^iv < 0.98 fm , 0.19 fm < Im a^^v < 0.37 fm . (27) 

In most recent publications, the value used for Im o^at is around 0.3 fm. However, for 
Rea^AT the estimates are still very different (compare, for example, Refs. |^ and [0]). In 
the present work we assume that 

a^^ = (0.55 + i0.30) fm . (28) 



The T-matrix t^'^ constructed in this way, reproduces the scattering length ( pSj) and the Su 
pole, but apparently violates the two-body unitarity since it does not obey the two-body 
Lippmann-Schwinger equation. 



2. Version II 



An alternative way of constructing the two-body T-matrix t"^^ is to solve the correspond- 
ing Lippmann-Schwinger equation with an appropriate separable potential having the same 
form-factors (^). However, a one-term separable T-matrix obtained in this way, does not 
have a pole at z = Eq — iT/2. To recover the resonance behavior in this case, we use the 
trick suggested in Ref. 0], namely, we use an energy- dependent strength of the potential 



Vik,k';z)=gik) 



A + C 



c 



9ik') 



where A is complex while C and ( are real constants. With this ansatz for the potential, 
the Lippmann-Schwinger equation gives the T-matrix in the form (^) with 



t{z) 



Ana"- 



A(c -z) + c(: 



I^vN J C~ z- [A(C - z)+ CC] /(I - i^J^I^/af 



(29) 



The constants A, C, and C can be chosen in such a way that the corresponding scattering 
amplitude reproduces the scattering length o^at and has a pole aX z = Eq — iV /2. 

This version of f'' also reproduces the scattering length ( pH]) and the Su pole. Moreover, 
it is consistent with the condition of the two-body unitarity. 
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3. Version III 

We can also construct the t^^ having the same form as version I, namely (^), with the 
same t{z) as in p^) but obeying the unitarity condition 

{l-27Tit'"'){l-2mt'T'y = 1 . (30) 

Of course, with the simple form (p3|), we cannot satisfy the condition (|30D at all energies. 
To simplify the derivations, we impose this condition on t^^ at z = Eq. Since Eq. (^) is 
real, it can fix only one parameter and we need one more condition to fix both the real and 
imaginary parts of the complex A. As the second equation, we used the imaginary part of 
Eq. ( p6D with a^iv given by (pSD. 

This procedure guaranties the two-body unitarity and gives the correct position of the 
resonance pole, but the resulting t^"^ gives ar,N which, of course, is different from the value 
(P5|), namely, it gives 

ttriN = (0.76 + i0.61) fm . (31) 

In what follows we use the three versions of the matrix t"^^ described above. All of them 
have the same separable form (^) but different t{z). A comparison of the results obtained 
with these three T-matrices can give us an indication of the importance of the two-body 
unitarity in the photoproduction processes. 



B. Photoabsorption N{r], 
In constructing the photoabsorption T-matrix t"'^, the Su dominance in the near- 



threshold region also plays an important role. It was experimentally shown that, at 



low energies, the reaction (0) proceeds mainly via formation of the S'n-resonance, i.e. 

-f + N—^Sn^N + T]. (32) 

This implies that in this energy region can be written in a separable form similarly to 
(p3|). To construct such a separable T-matrix, we use the results of Ref. where the 
was considered as an element of a multi-channel T-matrix which simultaneously describes 
experimental data for the processes 

n + N TT + N, Tc + N 1] + N, 
J + N IT + N, 'J + N r] + N 

on the energy shell in the S'n-channel. In the present work, we take the T-matrix t1^{E) 



from Ref. ||3^ and extend it off the energy shell via 



E) = ^;±|! t:Z(E) , (33) 

where k is a parameter. The Yamaguchi form-factors used in this ansatz, go to unity on 
the energy shell. Since k is not known, this parameter is varied in our calculations within a 
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reasonable interval 1 fm~^ < k < 10 fm~^ which is a typical range for meson-nucleon forces. 
It is known that f^^ is different for neutron and proton. In this work we assume that they 
have the same functional form and differ by a constant factor, 



Multipole analysis [0 gives for this factor the estimate A = —0.84 ± 0.15 . 

C. Nuclear subsystem 

Since the T-matrices and f'^ are poorly known and their uncertainties significantly 
limit the overall accuracy of the theory, it is not necessary to use any sophisticated ("real- 
istic") potential to describe the A^A^ interaction. Therefore we may safely assume that the 
nucleons interact with each other only in the S'-wave state. 

To obtain the necessary nuclear wave function ipQ, we solve the few-body equations of 
the Integro-Differential Equation Approach (IDEA) |^8|,^ with the Malfiiet-Tjon poten- 



tial . This approach is based on the Hyperspherical Harmonic expansion method applied 
to Faddeev-type equations. In fact, in the case of 5- wave potentials, the IDEA is fully equiv- 
alent to the exact Faddeev equations. Therefore, the bound states used in our calculations, 
are derived, to all practical purposes, via an exact formalism. 



IV. RESULTS AND DISCUSSION 

Figures 1-5 show the results of our calculations for the total cross section 

a = 

of the coherent process (|I]). The calculations were done for two nuclear targets, and ^He, 
using the three versions of described in the previous section. The curves corresponding 
to these three T-matrices are denoted by (I), (II), and (III), respectively. 

As can be seen in Fig. 1, the two versions of t^^, (I) and (II), give significantly different 
results despite the fact that both of them reproduce the same o^at and the S'n-resonance. 
This indicates that the scattering of the rj meson on the nucleons (final state interaction) 
is very important in the description of the photoproduction process. This conclusion is 
further substantiated when our curves are compared to the corresponding points (triangles) 




calculated for the ^He target in Ref. |^ where the final state interaction was treated using 
an optical potential of the first order. It is well-known that the first-order optical theory 
is not adequate at the energies near resonances. This is the reason why the calculations of 



Ref. ^1|] underestimate a near the threshold where with a^^v = (0.55 + ^0.30) fm the systems 



and 77 ^He have resonances ||23|| . 



A significant differences between the corresponding curves (I) and (II) in Fig. 1 imply 
that two-body unitarity is important as well. To clarify this statement, we compare in Fig. 
2 the results corresponding to the three choices of t{z) in (^). Surprisingly, the curves (II) 
and (III) almost coincide despite the fact that they correspond to different a^Ar while both 
obey the two-body unitarity condition. 
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The last three figures, Figs. 3, 4 and 5, show the dependence of the results on the choices 
of the parameters k and A = fj^ /f^ . Since nothing is known about k, we assume k, = a as 
the basic value for it. This can be motivated by the fact that both the elastic scattering and 
photoabsorption (production) of the 77 meson on the nucleon go via formation of the same 
Sii resonance. This means that at least one vertex, namely, r]N <-> S\\ should be the same 
for both the clastic scattering and photoabsorption. To find out how crucial the choice of 
K is, we did two additional calculations with k = 1 fm~^ and k = 10fm~^ (sec Fig. 3). We 
see that even with this wide variation the curves are not far from each other especially in 
the immediate vicinity of the threshold energy. Therefore the dependence on k is not very 
strong and the choice k — a can give us a reasonable estimate for the photoproduction cross 
section. Fig. 3 also shows an interesting tendency of increasing a when the range of the 
interaction becomes smaller (when k, grows). 

Figs. 4 and 5 for the and '^He targets, respectively, show the dependence of cr on the 
choice of the parameter A. An interesting observation here is that the cross section for 77 
photoproduction is more sensitive to this parameter with the tritium rather than the ^He 
target. This means that between these two nuclei, the tritium is a preferable candidate for 
a possible experimental determination of the ratio A. 

The cusp exhibited by all the curves at the threshold of total nuclear break-up, refiects 
losses of the fiux into the non-coherent channel. In a sense, this is a refiection of the four- 
body unitarity which the FRA equations are consistent with. 
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FIGURES 
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FIG. 1. Cross section of the coherent ??-photoproduction on the and ^He targets, calculated 
with the two versions of which are denoted as (I) and (II) respectively. All curves correspond to 
fljjAT = (0.55 + i0.30)fm, k = a = 3.316 fm~^, and A = —0.84. The triangles represent the points 
calculated in Ref. [^T[ for the '^He target within the optical model. 
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FIG. 2. Cross section of the coherent ?7-photoproduction on '^He, calculated with the three 
versions of t^'^ which are denoted as (I), (II), and (III) respectively. All three curves correspond to 
K = a = 3.316 fm~^ and A = —0.84 . For the curves (I) and (II) the o^at is given by Eq. (p8[) while 
for the third curve by Eq. (31) 
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FIG. 3. Cross section of the coherent 77-photoproduction on '^He, calculated with the version 
(II) of f^"^ with three values of the parameter k. All three curves correspond lo A = —0.84 . 
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FIG. 4. Cross section of the coherent //-photoproduction on ^He, calculated with the ver- 
sion (II) of t^"^ with three different values of the parameter A. All three curves correspond to 
K, = a = 3.316 fm~^ . 
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FIG. 5. Cross section of the coherent 77-photoproduction on ^H, calculated with the ver- 
sion (II) of t^"^ with three different values of the parameter A. All three curves correspond to 
K, = a = 3.316 fm~^ . 
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